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Abstract — Opportunistic scheduling (OS) schemes have been 
proposed previously by the authors for multiuser MIMO-SDMA 
downlink systems with linear combining. In particular, it has been 
demonstrated that significant performance improvement can 
be achieved by incorporating low-complexity linear combining 
techniques into the design of OS schemes for MIMO-SDMA. 
However, this previous analysis was performed based on the effec- 
tive signal-to-interference ratio (SIR), assuming an interference- 
limited scenario, which is typically a valid assumption in SDMA- 
based systems. It was shown that the limiting distribution 
of the effective SIR is of the Frechet type. Surprisingly, the 
corresponding scaling laws were found to follow elogA' with 
< e < 1, rather than the conventional log log K form. 

Inspired by this difference between the scaling law forms, 
in this paper a systematic approach is developed to derive 
asymptotic throughput and scaling laws based on signal-to- 
interference-noise ratio (SINR) by utilizing extreme value theory. 
The convergence of the limiting distribution of the effective SINR 
to the Gumbel type is established. The resulting scaling law 
is found to be governed by the conventional log log K form. 
These novel results are validated by simulation results. The 
comparison of SIR and SINR-based analysis suggests that the 
SIR-based analysis is more computationally efficient for SDMA- 
based systems and it captures the asymptotic system performance 
with higher fidelity. 

I. Introduction 

Opportunistic scheduling (OS) has recently attracted consid- 
erable research interest as a promising technique to improve 
system throughput by exploiting multi-user diversity with 
limited channel feedback [10]. Generally speaking, existing 
OS schemes can be classified into two categories, namely 
the time-sharing (TS) [10] and space-division multiple access- 
based (SDMA)-based [9] OS schemes. In TS-OS, only the 
mobile terminal (MT) with the best instantaneous channel 
conditions is scheduled in one slot, regardless of the number of 
beams employed by the base station (BS). In contrast, SDMA- 
based OS serves multiple MTs simultaneously with multiple 
orthonormal beams in each slot. Denote by M and N the 
number of transmit and receive antennas, respectively. It has 
been shown recently that the sum-rate of SDMA-based OS 
grows linearly with M whereas that of TS-OS increases only 
linearly with min(M, N) [8]. In addition to the more rapidly 
growing scaling law, SDMA-based OS is particularly attractive 
for practical systems with stringent latency requirements. 
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The SDMA-based OS in [9] was originally developed for 
systems with single-antenna MTs. For MTs with multiple 
receive antennas, [9] proposes to let each antenna compete for 
its desired beam as if it were an individual MT. As a result, 
each beam is assigned to a specific receive antenna of a chosen 
MT. Since signals received from the undesignated antennas of 
a chosen MT are discarded, this leads to inefficient utilization 
of multiple receive antennas. In [4], various linear combining 
techniques exploiting signals received by all receive antennas 
were proposed. The enhanced effective SINR is employed as a 
scheduling metric. Both analytical and simulation results in [4] 
have demonstrated that the system sum-rate performance can 
be significantly improved by using such combining techniques. 
For instance, the optimal combining technique can provide 
over 40% sum-rate improvement compared to the selection 
combining technique for M — 4 and N = 2 [4]. 

The theoretical analysis in [4] has been conducted based 
on SIR, assuming an interference-limited environment. The 
resulting scaling laws have a distinctive form, i.e. elogA" with 
< e < 1, which is very different from the conventional 
form log log AT derived based on signal-to-noise-ratio (SNR) 
[10] or SINR [8], [9] in the literature. Similar results have 
been independently developed for multicell systems in [2]. In 
this work, we introduce a systematic approach for deriving 
asymptotic throughput and scaling laws using SINR. The 
proposed approach stems from extreme value theory [3]. We 
prove that the cumulative distribution functions (CDFs) of the 
effective SINR obtained with linear combining converge to the 
Gumbel-type limiting distribution. Furthermore, we show that 
the SINR-based scaling laws for the proposed opportunistic 
beamforming and scheduling schemes follow the conventional 
log log A' form. Through comparison between the SIR and 
SINR-based analysis, it is argued that the SIR-based analysis 
is more computationally efficient for SDMA-based systems, 
and subsequently more effective in capturing the high-order 
behavior of the asymptotic system performance. To make 
comparison with our previous SIR-based analysis reported in 
[4], we concentrate on a practical system with M = 4 and 
N = 2 in this work. However, it should be emphasized that the 
analysis can be easily generalized for systems with arbitrary 
M and N. 

Notation : Vectors and matrices are denoted by boldface 
letters. ||-|| represents the Euclidean norm of the enclosed 
vector and | • | denotes the amplitude of the enclosed complex- 
valued quantity. Jjy is the N x N identity matrix. We use 
E {•} for expectation. Finally, log and In are the logarithms 
to the base 2 and e, respectively. 



II. Signal Model 



III. OS with Linear Combining 




Fig. 1 . A block diagram of the opportunistic MIMO SDMA downlink system 
under consideration. 



We consider the opportunistic MIMO-SDMA downlink 
system depicted in Fig. Q] where the BS is equipped with 
M transmit antennas and each of the K MTs has N receive 
antennas with N < M. Let {a m ; m = 1, 2, • • • , M} be a 
vector set containing M orthornormal beamforming vectors 
of length M . We focus on a particular time slot during which 
a beamforming vector set {a m } has been chosen from a 
common codebook shared by the BS and MTs. During the 
p-th slot, the transmitted signal can be expressed as 



x{p) 



M 

E 



= As(p), 



(1) 



where A = foi, 02,- • • , clm] is the unitary beamforming ma- 
I M and s(p) = [si(p), s 2 (p), s M (p)] T 
1 is the data vector transmitted in the 
p-th slot. The corresponding received signal by the fc-th MT 
can be written as 



[ax, a 2 , ■ 
trix with A H A 

with E{\s m (p)\ 2 } 



Vkip) = \/PkH k (p)x(p) + n k (p), 



(2) 



where H k is the channel gain matrix between the BS and 
the fc-th MT with independent and identically-distributed 
(i.i.d.) Rayleigh-distributed complex entries. Furthermore, the 
noise term n k (p) is modeled as CAf(0,I]y) and pf. is 
a constant related to the average received SNR given by 
E [ Pk \\H k {p)x{p)\\ 2 ) = Pk M. 

To keep our following analysis tractable, we concentrate 
on a homogenous system with p k = p in this work. For 
notational simplicity, we drop the temporal index p in the 
sequel. Furthermore, we refer to the SINR obtained by linearly 
combining signals from all receive antennas as the effective 
SINR in order to distinguish it from the observed SINR 
without combining. 



In this section, we briefly review the beamforming and 
scheduling schemes for MIMO-SDMA systems with linear 
combining techniques proposed in [4]. As shown in Fig.Q] in 
the beginning of a time slot, each MT evaluates the effective 
SINR for each beam by linearly combining the received signals 
with one of the following three combining techniques, namely 
selection combining (SC), maximum ratio combining (MRC) 
and optimum combining (OC) before returning the information 
about M effective SINRs to the BS. Note that OC performs 
active interference suppression by exploiting the interference 
structure, whereas MRC and SC simply intend to amplify the 
desired signal. It will be shown later that this characteristic 
interference-suppression feature of OC enables the scheduling 
scheme with OC to considerably outperform those with SC 
and MRC. 

Upon receiving the effective SINR information from all 
MTs, the BS schedules and starts data transmission to multiple 
MTs with the largest effective SINRs on different beams 
until the end of the current time slot. At each chosen MT, 
received signals from all antennas are linearly combined using 
one of the above linear combining techniques, followed by 
data detection. It is worth noting that the probability of 
awarding multiple beams to the same MT is rather small, as the 
number of MTs is large. Furthermore, recall that the minimum 
mean squared error (MMSE) and zero-forcing (ZF) receiver 
structures for MIMO receivers amount to combiners using OC 
and MRC for each beam, respectively. As a result, for an MT 
assigned with multiple beams, it can focus on one assigned 
beam at a time using the chosen combining technique while 
regarding all other beams as interfering sources. 

IV. SINR Analysis 

Define 7^ = max (7i, m , 72,™, • ■ • ,lK,m), for m = 
1,2, ■ ■ ■ , M. Assuming ^ ktm for fc = 1,2, - ■ ■ ,K, are i.i.d. 
with CDF Fx{x), the resulting average system throughput can 
be computed as [9]: 



c = iJ]T log (i + 7r *„)l 

lm=l J 



Ml log (1 -Fx) d[F x {x) 
10 



K 



(3) 



In the following, we first derive Fx(x) based on dif- 
ferent linear combining techniques before establishing their 
corresponding limiting distributions, i.e. lim [Fx(%j\ ■ By 

K— >oo 

exploiting the limiting distributions, we derive the asymptotic 
throughput and the corresponding scaling laws. In the sequel, 
we focus on a practical system with M = 4 and N = 2. 
However, it has been shown in [5] that the analysis can be 
easily generalized for systems with arbitrary M and N. 

A. Selection Combining (SC) 

We begin with the selection combining. Denote by x the 
maximum of the two SINR values of the z-th beam perceived 



by the two antennas at the k-th MT. The CDF of x can be 
derived based on the results in [9] and reads 

2 



(SC)/ s _ 



-x/p 



1 - 



(1 + x) 
respect 

the corresponding probability density function (PDF) 



(4) 



Differentiating F x c \x) with respect to x, we can obtain 



(SC)/ s _ 



1 - 



,-x/p 



(1 



(1 + x) 4 ' 



(5) 



It is straightforward to show that F x (x) and f x (x) 
satisfy the following equation 



/r(-) 



lim 



= P>0, 



(6) 



which is the necessary and sufficient condition for the lim- 



iting distribution of 



being of the Gumbel type 



[3]. Consequently, Fx {K) {x) = [Fx{x)] K converges to the 
following Gumbel-type distribution [3]. 



or equivalently, 



USC) 



? (SC) 

X (K) 
-,< sc > c n A h> sc > 



F^_[. (x) 



x > 



x > 0, 



(7) 



(8) 



where a { K and 6^? are normalizing factors affecting the 
shape and location of the limiting distribution, respectively. 
From extreme value theory, 6 ( ^ C) can be computed from the 
characteristic extreme of (0]l as [3] 

l-Ff C) (iC) = l (9) 



K 



Since (O is an exponential-linear equation of it is non- 
trivial to obtain the exact solution of b'^ ] in closed form. 

(SO 

Fortunately, since 1 — F x £ monotonically decreases from 
1 to whereas l/K € [1,0) for K = 1,2, • • • , oo, there 
always exists a unique solution of (O. Thus, we can resort to 
numerical methods to compute the numerical solution of 
It should be emphasized that 1 - F^ C) (b'^') tends to 1 as K 
approaches infinity, which implies that bf^ ^ increases with K. 

Furthermore, a'J5' can be obtained from solving the follow- 
ing equation. 



„(SQ _ TTl(SC) 

U K — e X 



(10) 



Similar to & ( JP, we can show that there always exists a unique 
solution of a ( J?. Therefore, the numerical solution of a'jJ 1 can 
be found by resorting to numerical methods. 

Finally, the throughput obtained with SC can be computed 
by substituting ([8]l into <[3j and reads 



MSC) _ _t_ 
In 2 



1 - e 



1 + x 



■ dx. 



(11) 



Let z = e 



and £ = exp(6f /a iS £). We have x 



— a'^'lnz and dx — dz. Thus, ( fTTT i can be rewritten as 



4 
ln~2 

4 
W2 



1 



K 



1 - 

4 

„ (i 



af> lnz 



(12) 
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a'^lnz) z 



dz- 



4 (1 



A' 



In z) z 



(13) 



The limit of the first term on the right-hand-side (R.H.S) 

4 

of (fL3l becomes negligibly small as lim 7 = while the 

A— >oo £ 

limit of the second term can be computed by exploiting the 
approximation of a c J? ~ p as follows. 



lim 



'A 



dz 



lim 4 log (6f 



K-»oom2,/4 (1 - a^'lnz) z 
Thus, the corresponding scaling law is given by 



lim 



(7(SC) 



a— 41og(&f>) 



1. 



(14) 



(15) 



In particular, for p = 1, we can approximate bg? and a'J? as 
ln2if -21n(l + ln2if), (16) 



L(SC) 



'A 



1, 



(17) 



respectively. 

Subsequently, the scaling law can be written as follows. 



lim 



= 1, 



(18) 



a^oo 4 log (In 2K - 2 In (1 + In 2if )) 
which stands for a typical scaling law in the loglogi^T form. 

B. Maximum Ratio Combining (MRC) 

The effective SINR obtained with MRC can be expressed as 
a ratio of two random variables given by x = Yfp+y 7 wnere 
z and y are \ 2 distributed random variables with 2N and 
2M—2 degrees of freedom corresponding to the instantaneous 
signal power of the desired signal and the interfering signal, 
respectively. In particular, for M = 4 and N = 2, we have 
[6], [7] 



-x/p 



Qxe- X ' p I2xe~ x / p 



/•(MRC) / \ _ 

Jx { > ~ p 2 {l + xf p{l + xY ' (1 + a;) 5 
and the corresponding CDF can be expressed as 

e -x/ P xe - x /P 3xe~ x /P 



F^ix) 



1 



(1 + x) 
It can be shown that 
1 



lim 



p(MRC) 



(X) 



+ x y (l + xY 



= p > 0. 



(19) 



(20) 



(21) 



Therefore, the limiting distribution of 
of the Gumbel type. 



F^ix) 



A 



is also 



Following similar steps as in the previous section, we have 

,(MRC) 



(MRC! 1 (MRC) 



and 



/^(MRC) _ ^ 

~ In 2 



lim 



l-e 



■ dx. 



£<(MRC) 



^oo41og(6^) 



= 1, 



(22) 



(23) 



where a'^ RCI and 6 ( ™ c) are the corresponding normalizing 
factors. In particular for p = 1, we can show that 



a' 



K 

(MRC) 

K 



In 3K - 2 In (1 + In K) , (24) 
1 (25) 



(26) 



and the scaling law has the following log log K form. 

p(MRC) 

k™L 4 log (In 3K - 2 In (1 + In K)) 
C. Optimal Combining ( OC) 

The CDF of the effective SINR obtained OC using N 
receive antennas in the presence of M — 1 interfering sources 
has been derived in [1]. For M = 4 and N — 2, the 
corresponding CDF takes the following form. 



c >(*) = 1 



,-x/p 



{l + xf (1 
and the corresponding PDF is 

-x/p 



x) 3 p(l + x) 3 



(27) 



/? C) (*) 



xe 



Since lim 



p 2 (l + x) 
1 - Ff*(x) 



™ jf C) (z) 



j [(3p + 1) x + (6p 2 + 6p + 1)] . 

(28) 

p > 0, the limiting distribution 



of 



F^(x) 



K 



is also of the Gumbel type. Similar to the 



cases of SC and MRC, we can show that 



(OC) T (OC) 



r (OC) _ * 



■ dx. 



and 



C(oo 



*r lnn oo41og(&'° C| ) 



1. 



(29) 



(30) 



where and 6^ c) are the corresponding normalizing factors. 
In particular for p = 1, we can show that 



UOCj 



'K 



In IK- 2 In In K, 
1 



(31) 
(32) 



and the sum-rate scales like the following log log K form. 



r<{oc) 
<- P =i 



lim p — = 1. 

k^oo 4 log (In AK - 2 In In K) 



(33) 




Fig. 2. Simulation versus analytical results with numerical normalizing 
factors for p = 1, 5. 



V. Simulation Results 

In this section, simulation is performed to confirm our SINR 
analysis derived in Sec. [IV] Unless otherwise specified, we set 
M = 4 and N = 2. 

We first compare the asymptotic throughput shown in ( fTTT l. 
(l22i > and ( |29l against their corresponding simulation results. 
Figure [2] shows the asymptotic throughput curves using the 
numerical normalizing factors obtained by numerical methods 
for p = 1 and 5. Inspection of Fig.|2]reveals that the analytical 
results shown in (fTTI i. (l22l and (1291 are in accord with the 
simulation results. Despite that the asymptotic analysis is 
achieved by assuming a large K, Fig. [2] indicates that the 
asymptotic analysis is also very accurate for smaller K values. 
Furthermore, Fig. [2] confirms that the scheduling scheme with 
OC can substantially outperform those with MRC and SC 
whereas the improvement provided by MRC is more apparent 
in the presence of stronger noise. This is because the scheme 
with OC is designed to maximize SINR whereas MRC intends 
to maximize SNR. 
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Fig. 3. Simulation versus analytical results with approximated normalizing 
factors for p = 1. 



Next, rather than the numerical solutions, Fig. [3] depicts the 



average sum-rates using the approximated normalizing factors 
computed in ( TToT ). ( f24T > and (l3TT l together with a,K ~ 1 for 
p = 1. Since the approximation expressions have been derived 
by assuming a large K, the analytical curves shown in Fig. [3] 
approach the simulated curves only when K becomes large. 
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Fig. 4. Comparison of the numerical and approximated solutions of the 
normalizing factors for p 1. 

Finally, to inspect the approximation accuracy of (TToT l. (|24| | 
and (|3TT i. Fig. 0] shows the numerical and approximated nor- 
malizing factors as a function of the number of MTs, K. Since 
solving the exact solutions to the normalizing factors involves 
the linear-exponential functions, it is in general non-trivial to 
obtain accurate closed-form expressions for the normalizing 
factors, which compromises the accuracy of the subsequently 
derived scaling laws. 

VI. Comparison between SIR and SINR analysis 

It is interesting to compare the SINR analysis derived in 
this work with our previous SIR analysis reported in [4]. 

1. ) On the one hand, it is easy to verify that the CDFs 
of the effective SINR in ©, d20ji and d221> converge the 
corresponding CDFs of the effective SIR reported in [4] as p 
tends to infinity, respectively. On the other hand, our SIR and 
SINR-based analysis suggests that the limiting distributions 
of the effective SIR and SINR do not belong to the same 
domain of attraction. Instead, they are of the Frechet-type and 
Gumbel-type, respectively. It is natural to conjecture that the 
limiting distribution function of SINR might also converge to 
the Frechet-type if the noise power becomes zero. However, 
our results reveal that this intuition is not true. This is because 
that the limit operator is not commutative in general. 

2. ) It is generally more difficult to obtain the normalizing fac- 
tors in the SINR analysis than the SIR analysis since the SINR- 
based analysis involves exponential-type CDFs and requires 
solving exponential-linear equations such as ©. Therefore, it 
is more computationally advantageous to derive the scaling 
laws in the SIR-based analysis compared to the SINR-based 
analysis in the presence of strong interference. 

3. ) When computing the normalizing factors in the SINR- 
based analysis, we have to carefully take into account the 



high-order terms in Fx(x). For instance, if the high-order 
terms in Fx (x) in (@J, (f20b and (l27l i are ignored, the resulting 
simplified CDFs for different schemes will all lead to the 
same set of normalizing factors, i.e f. , f ^ = 4r. Thus, 
the performance of OS schemes with different combining 
techniques cannot be distinguished based on their scaling laws. 
Since it is generally much easier to compute the normalizing 
factors with high accuracy in the SIR-based analysis [4], we 
argue that the SIR-based scaling laws can better characterize 
the actual performance of different OS schemes by focusing 
on the interference-limited scenarios. 

VII. Conclusion 

In this paper, we have developed a systematic approach to 
derive the SINR-based asymptotic throughput and scaling laws 
for OS schemes by utilizing extreme value theory. In particular, 
we have investigated the asymptotic throughput and scaling 
laws of the OS schemes proposed for MLMO-SDMA systems 
with different linear combining techniques. Our analytical 
results have shown that the limiting distribution of the effective 
SINR is of the Gumbel type and the scaling laws follow 
the log log K form. Simulation results have confirmed the 
effectiveness in improving system throughput by incorporating 
low-complexity linear combining techniques in OS schemes. 
Finally, based on the comparison of SIR-based and SINR- 
based analysis, we have argued that the SIR-based analysis is 
more advantageous in providing insights into the scheduling 
performance for SDMA-based systems. 
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